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J /ifj suppression at SPS and RHIC in the comovers approach 
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The NA50 collaboration data on the J/il) suppression are compared with the results 
obtained in a comovers approach based on the Dual Parton Model (DPM). Predictions 
for the J / ip suppression versus the charged multiplicity - measured in the rapidity region 
of the dimuon trigger - are given for SPS and RHIC energies. 

The NA50 collaboration has observed an anomalous J/ip suppression in Pb Pb colli- 
sions i.e. a suppression that exceeds the one expected from the extrapolation of the pA 
and SA data, and exhibits an interesting centrality pattern. These data [|l| have been 
interpreted either as the result of one (or two) deconfining phase transition(s) or as due 
to the interaction with comovers. Our aim is to examine to what extent the observed 
pattern can be reproduced in the latter approach. Here, the J/^p survival probability is 
the product of two factors Sabs{b,s) ■ Sco{b,s). The first factor represents the suppres- 
sion due to nuclear absorption of the cc pair. Its expression, given by the probabilistic 
Glauber model, is well known. It contains a parameter, the absorptive cross-section aabs- 
The second factor Sco{b, s) represents the suppression resulting from the interaction with 
comovers. Its expression [0-H] is 



Sco{b, s) = exp <^ -aco Ny^^{b, s) In 
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Here iV™ represents the density of comovers (positive negative and neutrals) in the rapidity 
region of the dimuon trigger {0 < y* < 1). We determine it in the Dual Parton Model (see 
Eq. (^) below). Nf = 1.15 fm~^ is the corresponding density in pp, and dco is an averaged 
interaction cross-section. In the calculations below we take aco = 1 mb and aabs = 4.5 mb. 

In this way we can compute the J/ip suppression at each impact parameter. However, 
experimentally, the ratio of J/i/) over DY is plotted as a function of either Et or the energy 
of the zero degree calorimeter Ezdc- Et is the transverse energy of neutrals deposited in 
the NA50 calorimeter, located in the backward hemisphere (1.1 < yiab < 2.3). Using the 
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proportionality between Et and multiplicity, we have 

ETib) = lqN-Jb) . (2) 

Here the multiplicity of comovers is determined in the rapidity region of the NA50 
calorimeter. The factor 1/3 is introduced because only the energy of neutrals is recorded. 
In this way q is close to the transverse energy per particle, but it contains also the cal- 
ibration of the calorimeter. From a fit to the Et distribution we obtain q = 0.6 GeV. 
An alternative determination of q is obtained from the measured correlation between Et 
and EzDC (see below). The latter is defined as EzDcip) = [A — n^(6)]i5j„. Here ua 
is the number of participants of nucleus, A — ua the number of spectators, and Ein = 
158 GeV is the beam energy. We see from Eqs. (|l|) and @ that, in order to describe the 
centrality dependence of the J/ ip suppression, it is paramount to have a good description 
of the h dependence of Ny" - both in the rapidity region of the dimuon trigger and in the 
one of the Et calorimeter. The 6-dependence of charged multiplicity at mid-rapidities 
has been measured at SPS and RHIC. It has been shown in that these data are well 
described in DPM, where the multiplicity is given by a linear combination of the number 
of participants and the number of binary collisions 

N;\b) = Ayib) UAib) + Byib) nib) (3) 

with coefficients that depend on h and y and can be calculated in the model. Their 
values at mid-rapidities are given in 0. Using Eq. (j^), with coefficients calculated in 
the rapidity region of the calorimeter, we obtain a good description of the measured 
Et — EzDC correlation with q = 0.6 GeV (Fig. 1)Q The correlation Et — & is given by 
P{Et, b) oc exp{-{ET - ET{b)f /2qaET{b)} with q = 0.6 GeV and a = 0.88. 

The model allows to compute the ratio J/tp over DY versus either Et or Ezdc from 
peripheral collisions up to the knee of the Et distribution. To go beyond it, we have to 
introduce the fluctuations responsible for the tail of the Et distribution (Eq. @) 
gives only the average value of Et at each b) . They have been introduced in the model 
by multiplying Ny°^^{b) in Eq. (||) by F{b) = ET/ET{b), where Et is the measured value 
of the transverse energy. 

The results are presented in Fig. 2a and compared with the NA50 data. We see that 
from peripheral collisions up to the knee of the Et distribution, the data are well described. 
However, beyond the inflexion point at the knee, the decrease in the data is sharper than in 
the model. Note, however, that the data beyond the knee are obtained with the so-called 
minimum bias (MB) analysis. Only the ratio J/ip over MB is measured and it is multiplied 
by a theoretical ratio DY/MB. In the model, this ratio (as well as the ratio J/ip/MB) is 
essentially flat beyond the knee - due to the fact that the tail of the Et distribution of 
hard {J/ip, DY) and soft (MB) processes is assumed to be the same. In sharp contrast 

"'^In 1^-^ an approximation was made in which the dependence on b of the coefficients A and B was 
dropped. Moreover, a two-string approximation was used for each A^iV-coUision. It turns out that these 
approximations have a very small effect on the ratio J/iJj over DY at each b. However, they do change 
the correlation between b and Et - or between Et and Ezdc- As a consequence, an extra term in the 
r.h.s. of Eq. (||) was needed in |^,^ in order to reproduce this correlation. 
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with this behaviour, the measured ratio J/ip over MB has a very pronounced turn-over 
at the knee. It is most important to determine whether such a turn-over is due to an 
increase in the J/ip suppression, or, on the contrary, to a small difference in the widths 
of the tails of the distributions of soft and hard processes. In the latter case a similar 
turn-over would also be present in the ratio DY/MB - resulting in a flatter behaviour of 
the ratio J/ip over DY beyond the knee. On the theoretical side the calculation beyond 
the knee is also subject to uncertainties. In particular, the rapidity regions of the dimuon 
trigger and the Ej- calorimeter are far apart and the fluctuations in the two regions could 
be different. In view of that, it would be very interesting to measure the J/-?/' suppression 
versus charged multiplicity - with the dimuon and multiplicity triggers sitting in the same 
rapidity region. The results of the model for SPS and RHIC are shown in Fig. 3. 

The NA50 collaboration has presented new data on the J/ip suppression versus Ezdc- 
In so far as a theoretical model describes well the Et — Ezdc correlation, a set of data can 
be plot versus either variable without further constraints on the model. However, the new 
data are important since they include many peripheral points. The results of the model 
are presented in Fig. 2b. The centrality dependence of the J/ip suppression is reasonably 
well reproduced. However, the absolute normalization used in Fig. 2a is here 10 % too 
large. (Note that the normalization of the data in Fig. 2b is not measured. It has been 
fixed [m from the one in Fig. 2a using only points in the range 60 < Et < 100 GeV). This 
point has important consequences for the interpretation of the data. Indeed, it is obvious 
from a comparison of figures 2a and 2b that, if the relative normalization of the two sets 
of data were correct, the onset of the anomalous suppression in Ezdc would take place 
at a value of b significantly larger than in E^- This problem is solved with the change of 
normalization discussed above. However, in this case there would be several experimental 
points sitting above the NA50 nuclear absorption curve. More important, independently 
of their absolute normalization, the data of Fig. 2b for Ezdc < 28 TeV exhibit a. J/ip 
suppression significantly steeper than the NA50 absorption model - indicating that the 
anomalous suppression is already present in very peripheral collisions, where the density 
is significantly lower than the maximal one in S-U. 

In conclusion, the NA50 data on the J/ip suppression from peripheral collisions up to 
the knee of the Et distribution can be described in a comovers approach. (In contrast, 
both the NA50 absorption model and the model fail to reproduce the peripheral data). 
Beyond the knee the obtained suppression is too small. However, in this region the 
theoretical uncertainties are large and the data contain a theoretical input that has to be 
checked experimentally. 
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Figure 1. Ex — Ezdc correlation. Figure 3. Ratio J /tl) over DY normalized 

to pp versus charged multiplicity in Pb- 
Pb collisions at piab = 158 AGeV/c, both 
in the range < < 1, (full line) and 
at ^ = 200 AGeV in the range -0.5 < 
y* < 0.5 (dotted line). 
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Figure 2. a): Ratio J/ip over DY versus Et compared to NA50 data The full curve is 
the theoretical prediction, b): Ratio J/ip over DY versus Ezdc compared to preliminary 
results presented by NA50 The full line is obtained from the full line in a) using the 
calculated Et — Ezdc correlation and changing the normalization by a factor 0.92 (see 
text). The dashed line is obtained computing for each 6, the value of Ezdc and the value 
of the ratio R - with the same change in normalization. In both figures the dotted line is 
the NA50 absorption model, fitting pA and SU. 



